Direct microscopic calculations of nuclear level densities in the shell model Monte 

Carlo approach 
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Nuclear level densities are crucial for estimating statistical nuclear reaction rates. The shell model 
Monte Carlo method is a powerful approach for microscopic calculation of state densities in very 
large model spaces. However, these state densities include the spin degeneracy of each energy level, 
whereas experiments often measure level densities in which each level is counted just once. To enable 
the direct comparison of theory with experiments, we introduce a method to calculate directly the 
level density in the shell model Monte Carlo approach. The method employs a projection on the 
minimal absolute value of the magnetic quantum number. We apply the method to nuclei in the 
iron region as well as the strongly deformed rare-earth nucleus Dy. We find very good agreement 
with experimental data including level counting at low energies, charged particle spectra and Oslo 
method at intermediate energies, neutron and proton resonance data, and Ericson's fluctuation 
analysis at higher excitation energies. 

PACS numbers: 21.10.Ma, 21.60.Cs, 21.60.Ka, 21.60.De 



Introduction. The level density is among the most im- 
portant statistical properties of the atomic nucleus. It 
appears explicitly in Fermi's golden rule for transition 
rates and in the Hauser-Feshbach theory [1| of statis- 
tical nuclear reactions. Yet its microscopic calculation 
presents a major theoretical challenge. In particular, cor- 
relations have important effects on nuclear level densi- 
ties but are difficult to include quantitatively beyond the 
mean-field approximation. The configuration-interaction 
(CI) shell model is a suitable framework to include both 
shell effects and correlations. However, the dimension of 
the required model space increases combinatorially with 
the number of single-particle states and/or the number of 
nucleons, and conventional shell model calculations be- 
come intractable in medium-mass and heavy nuclei. This 
difficulty has been overcome using the shell model Monte 
Carlo (SMMC) approach d-Q. The SMMC has proved 
to be a powerful method to calculate microscopically nu- 
clear state densities [6l4ll|. 



The SMMC method is based on a thermodynamic ap- 
proach, in which observables such as the thermal energy 
are calculated by tracing over the complete many-particle 
Hilbert space. Thus, the SMMC state density takes into 
account the magnetic degeneracy of the nuclear levels so 
that each level of spin J is counted 2J + 1 times. 

However, experiments often measure the level density, 
in which each level is counted exactly once, irrespective of 
its spin degeneracy 12|, Il3l ] . To make direct comparison 
of theory with experiments, it is thus necessary to be able 
to calculate the level density within the SMMC approach. 
A spin-projection method, introduced in Ref. [lOl, can be 
used to calculate the level density pj{Ex) for a given 
spin J at each excitation energy E^- While the state 
density is given by p{Ex) — X]j(2>/ + 1)/5j(^k), the total 
level density is p{Ex) —J2j Pj{Ex)- However, this latter 



formula is not useful for practical calculations because 
the statistical errors of pj{Ex) increase with J, and the 
resulting statistical errors in p{Ex) are too large. 

Here we introduce a simple method to calculate di- 
rectly and accurately the level density in SMMC. We 
present level density calculations of medium-mass nuclei 
in the iron region and of the well-deformed nucleus -'^^■^Dy. 
We find good agreement with a variety of experimental 
data including level counting at low energies, charged 
particle spectra and Oslo method at intermediate ener- 
gies, neutron and proton resonance data, and Ericson's 
fiuctuation analysis at higher excitation energies. We 
note that our method can be applied more generally to 
many-particle systems with good total angular momen- 
tum. 

Level density in SMMC. We make the observation that 
for any nuclear level with spin J and magnetic quantum 
number degeneracy of 2 J -I- 1, the state with the lowest 
possible non-negative spin projection M appears exactly 
once. Denoting by pM the level density for a given value 
of the spin projection M, the total level density for even- 
even and odd-odd nuclei (whose spin is integer) is given 
by /5 = pm=Oi while for odd-even nuclei (whose spin is 
half-integer), the total level density is p = Pm=i/2- 

The M-p rojected level density can be calculated as 
in Ref. llOl For a nucleus described by a shell model 
Hamiltonian H and at inverse temperature /3 — 1/T, the 
SMMC method is based on the Hubbard-Stratonovich 
(HS) transformation jl4| e"^^ = / £> [ct] G^ C/^ , where G^ 
is a Gaussian weight and Ua is a one-body propagator de- 
scribing non-interacting nucleons in time-dependent aux- 
iliary fields (7. For a quantity X that depends on the 
auxiliary fields a, we define 



[D\a\W(a)X^ 



(1) 



where W{a) = Ga-\Tr Ua-\ is the weight used in the Monte 
Carlo samphng and ^^ = Tr C/cr/|Tr [7^1 is the Monte 
Carlo sign function. Here and in the following, the traces 
are evaluated in the canonical ensemble for fixed number 
of protons and neutrons, which in turn can be calculated 
from grand-canonical traces by particle-number projec- 
tion. 

The M-projected thermal energy Em {(3) = {H)ai is 
calculated using 
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The trace TtmX at fixed spin component M can be 
calculated by a discrete Fourier transform 
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where (pk (fc = — J^, . . . , J^) are quadrature points tpk 



'^l+TH ^^^ '^■5 ^^ ^^^ maximal spin in the many-particle 
shell model space. 

The M-projected canonical partition function Zm(/3) 
is calculated by integrating the thermodynamic relation 
-din Zm/cLP = EmW), taking ZmW = 0) to be the total 
number of levels with the magnetic quantum number M. 
For the lowest non- negative value of M, Zm{P — 0) is the 
total number of levels without counting their magnetic 
degeneracy. The M-projected level density pm{Ex) is 
then calculated in the saddle-point approximation 
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where Sm and Cm are, respectively, the M-projected 
canonical entropy and heat capacity 



Sm = In Zm + f^E^ 
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In the calculation of Cm we implemented the method of 
Ref. 1 151 . in which the numerical derivative is carried out 
inside the HS path integral. This enable us to take into 
account correlated errors, thus reducing significantly the 
statistical errors in the heat capacity compared to a direct 
numerical derivative of the thermal energy. Equation Q 
is analogous to the formula used for the state density [G] 
in which the corresponding quantities do not include M 
projection. 

The projection on the spin component M usually in- 
troduces a sign problem that leads to large fluctuations of 
observables at low temperatures (even for a good-sign in- 
teraction). However, for even-even nuclei Tr I e^'^'-'^-C/o- ) 

is almost always positive (for a good-sign interaction), 
and using Eq. Q with M = and X = f/o- we have 
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FIG. 1: Level densities versus excitation energy Ex for 
^''Fe, ""Ni, ^^^Ni and "^"Co. SMMC level densities p{Ex) = 
pM=o{Ex) (solid circles) are compared with various experi- 
mental data sets [I^l : level counting at low excitation energies 
(open diamonds), charged particle spectra [la ] at intermedi- 
ate energies (dashed lines), and Ericson's fluctuation analy- 
sis [iTJ] at higher energies (open circles). For ®°Co there is 
also the proton resonance data (open square) [la . [l9| . 



TrM=oC^(T > 0. Thus the level density of even-even nu- 
clei can be calculated accurately down to low excitation 
energies without a sign problem. 

Medium-mass nuclei. We demonstrate the SMMC cal- 
culation of level densities for medium-mass nuclei in the 
iron region using the CI shell model Hamiltonian of Ref. |6| 
in the complete pfgg/2 shell. 

In Fig. [T] we compare SMMC level density calculations 
(solid circles with error bars) for ^^Fe, ^"Ni, ^^Ni and 



^°Co with various experimental data compiled in Ref. 113: 
(i) level counting at low excitation energies (open dia- 
monds), (ii) charged particle reactions such as (a, a'), 
(p,p'), (p, a) and {a,p) at intermediate excitation en- 
ergies (dashed hnes) [l6J, and (iii) Ericson's fluctuation 
analysis at higher excitation energies (open circles) [l7[ . 
For ^°Co there is also high-resolution proton resonance 
data at around 8 MeV (open square) [l8|, |l9|. Overah, 
we find good agreement between the SMMC calculations 
and the experimental data. 
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FIG. 2; Top: The SMMC state density (open squares) and 
level density (solid circles) versus excitation energy Ex for 
®°Co. The experimental level density data foUows the same 
convention of Fig. [l] Bottom: thermodynamic moment of 
inertia for ^''Co extracted from the ratio of the state density 
to the level density (solid circles). The dashed line is the 
rigid-body moment of inertia. 



Spin-cutoff parameter. In the spin-cutoff model, the 
spin distribution pj{Ex) is given by 



Pj{E.,) 



'<->ai^»-^ 



(6) 



where p{Ex) is the total state density and Uc — adE^) is 
an energy-dependent spin-cutoff parameter. The distri- 
bution (|ni) is normalized such that J^ji'^'^'^ ^)pj{Ex) ~ 
p{Ex). Equation ([B]) can be derived in the random cou- 
pling model of individual spins [20|. In this model, the 
level density p{Ex) can be calculated to be 



p{Ex)=J2pjiE.) 



2TTa, 



-p{E.) 



(7) 



where the sum over spin is calculated by converting 
it to an integral. An effective spin-cutoff parame- 
ter can then be estimated from the ratio of the total 
state density to the total level density, i.e., adE^) = 
{2tt)~^^'^ p{Eri:) / p{Ex) ■ The spin-cutoff parameter can be 
converted to a thermodynamic moment of inertia I using 
al = IT/h\ 

We have extracted such moment of inertia / from the 
calculated SMMC state and level densities of ^^Fe and 
^°Co. In Figs. [2] and [3] we show the corresponding state 
densities (open squares) and level densities (solid circles) 
and the corresponding moment of inertia / (bottom pan- 
els) versus excitation energy Ex- For the odd-odd nu- 
cleus ^"Co the moment of inertia depends only weakly 



FIG. 3: As in Fig. H but for ^^^Fe. 



on excitation energy. However, for the even-even nucleus 
^^Fe we observe a suppression of the moment of inertia 
at low excitation energies. This reflects the reduction in 
the state-to-level density ratio that originates in pairing 
correlations, and is consistent with the results found in 
Ref. [lO| in which the moment of inertia was extracted 
from the spin distribution. 
Rare-earth nucleus ^^^Dy. 



In Refs. |21j and |2^ we ex- 
tended the SMMC approach to heavy nuclei in the rare- 
earth region using the 50-82 major shell plus the I/7/2 or- 
bital for protons, and the 82-126 major shell plus O/111/2 
and I59/2 orbitals for neutrons. We described success- 
fully the rotational character of the strongly deformed 
nucleus ^^^Dy 21] as well as the crossover from vibra- 
tional to rotational collectivity in families of samarium 
and neodymium isotopes [22| . 

Fig. H shows the SMMC level density p{Ex) = 
Pm=o{Ex) (soUd circles) and SMMC state density p[Ex) 
(open squares) for ^^^Dy. The SMMC level density com- 
pares well with various experimental data sets: (i) level 
counting (solid histograrns) 1231 . l2J|, (ii) renormalized 
Oslo data (open circles) 25L |26j and (iii) neutron reso- 
nance data (triangle) [27| . We find very good agreement 
between the various data sets and the SMMC level den- 
sity 

Conclusion. In conclusion, we have used a spin- 
component projection method to calculate directly and 
accurately the SMMC nuclear level density p{Ex) as the 
projected density pm=o{Ex) for even-even and odd-odd 
nuclei. The method is easily extended to odd-even nu- 
clei by using p{Ex) — Pm=i/2{Ex)- This method allows 
us to make direct comparison with experimental data. 
We find very good agreement between the SMMC level 
density and the experimental data for nuclei in the iron 
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FIG. 4: Level density and state density in ^'^^Dy. The SMMC 
level density (solid circles) is compared with the state den- 
sity (open squares). Also shown are experimental data sets 
level counting at low excitation ener- 
2JI, Oslo data at intermediate energies 
and the neutron resonance data (tri- 



for the level densit 
gies (histograrns) [2 




(open circles) 
angle) [13]. 



region and for the rare-earth nucleus ^^^Dy. 
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